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IC 310 has recently been identified as a gamma-ray emitter based on observations at GeV energies 
with Fermi-LAT and at very high energies (VHE, E > 100 GeV) with the MAGIC telescopes. 

Despite IC 310 having been classified as a radio galaxy with the jet observed at an angle > 10 
degrees, it exhibits a mixture of multiwavelength properties of a radio galaxy and a blazar, possibly 
making it a transitional object. On the night of 12/ 13 th of November 2012 the MAGIC telescopes 
observed a series of violent outbursts from the direction of IC 310 with flux-doubling time scales 
faster than 5 min and a peculiar spectrum spreading over 2 orders of magnitude. Such fast variability 
constrains the size of the emission region to be smaller than 20% of the gravitational radius of its 
central black hole, challenging the shock acceleration models, commonly used in explanation of 
gamma-ray radiation from active galaxies. Here we will show that this emission can be associated 
with pulsar-like particle acceleration by the electric field across a magnetospheric gap at the base 
of the jet. 


I. INTRODUCTION 


The nearby lenticular (SO, z = 0.0189) galaxy 
IC 310 located in the Perseus cluster exhibits an active 
galactic nucleus (AGN). This object has been detected 
at high energies (above 30 GeV) with Fermi /LAT [23| 
as well as at TeV energies IIQ. The jet of IC310, ex¬ 
tending in the outward direction from the center of 
the cluster led to early assi gnm ent of this object as 
a head-tail radio galaxy 19l 1271. [20l|. However, using 
the Very-Long-Baseline Interferometry (VLBI) tech¬ 
nique, a parsec-scale one-sided jet was found to follow 
the large-scale jet within about 10° [r|- The align¬ 
ment of the jet at different scales, without any hints of 


bending put in doubt the above classification. Instead, 
the inner jet appears to be blazar-like with a missing 
counter jet due to relativistically boosted emission. 
Further indications for transitional behavior between 
a radio galaxy and a blazar were found in IC310 in 
various energy ranges [2j|. The mass, M, of the black 
hole of IC310 can be inferred from its relation with 
the velocity dispersion, a, of the host galaxy [fl H5|. 
namely M = (3+2) x 10 8 Mq. 

MAGIC (Major Atmospheric Gamma Imaging 
Cherenkov) is a system of two 17-m diameter Imag¬ 
ing Atmospheric Cherenkov telescopes located on La 
Palma, Canary Islands. It allows observations of 
gamma-ray sources with energies above 50 GeV. Dur- 
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ing the observations of the Perseus cluster performed 
in the end of 2012 MAGIC telescopes revealed an ex¬ 
treme gamma-ray flare from IC 310 on the night of 
12/ 13 th of November [6|. In addition, the source was 
observed in radio band by European VLBI Network 
(EVN) during October/November 2012. 

In Section mi we report the data analysis and re¬ 
sults of the MAGIC observations during the flare and 
the radio observations. In Section HTH we discuss pos¬ 
sible interpretation of the ultrafast variability of the 
gamma-ray emission observed from IC310. 


II. RESULTS 

A. MAGIC 

MAGIC telescopes were observing the Perseus clus¬ 
ter on the night of 12/13 th of November for 3.7 h. 
The observations consisted of 4 pointings, two of the 
them with a standard offset of 0.4° with respect to 
IC 310 and the remaining ones are at a distance 0.94° 
away from the object. The signal extraction and cali¬ 
bration of the data, the image parametrization, the 
direction and energy reconstruction as well as the 
gamma-hadron separation were applied with the stan¬ 
dard analysis software MARS as described in [30| • 

In the night of the flare a strong signal of 507 
gamma-like events above 300 GeV in the region 
around IC 310 in excess of the background estimated 
as 47 events was observed. Due to still limited statis¬ 
tics of events and a very rapid variability behavior, the 
classical approach for the calculation of light curves 
in gamma-ray astronomy which is based on the fixed 
width of the time bins is not optimal in this case. We 
used instead a method similar to the one commonly 
used for data of X-ray observatories for the comput¬ 
ing of energy spectra. We first identify all periods in 
the data during which the telescopes were not opera¬ 
tional (in particular < 1 min gaps every 20 min when 
the telescope is slewing and reconfiguring for the next 
data run). Afterwards, we bin the remaining time 
periods based on a fixed number (in this case 9) of 
ON events per bin. We estimate the number of back¬ 
ground events in each time bin from four off-source 
regions at the same distance from the camera cen¬ 
ter. Using toy MC simulations we validated that this 
method limits the bias in flux value and its error .6j. 
As the signal to background ratio above 300 GeV is 
much larger than 1 this assures that the precision of 
individual points in the light curve is close to 3er. 

The resulting light curve is presented in Fig. [I] 
The mean flux above 300 GeV during this period is 
# mem = (6.1 ± 0.3) x 10 -11 cm _2 s _1 . This is four 
times higher than the high state flux of (1.60±0.17) x 
10 -11 cm _2 s _1 reported in [5|. The emission is highly 
variable, fitting the light curve in the full time range 



FIG. 1: Light curve of IC 310 observed with the MAGIC 
telescopes in the night of November 12/13 th , 2012, above 
300 GeV. As a flux reference, the two gray lines indicate 
levels of 1 and 5 times the flux level of the Crab Neb¬ 
ula, respectively. The precursor flare (MJD 56243.972- 
56243.994) has been fitted with a Gaussian distribution. 
The figure is reprinted from f6]. 



FIG. 2: Zoom of the first big flare seen in the light curve 
of IC 310 above 300 GeV. Black lines show exponential fits 
to the rising and decay edges to the substructures in the 
light curve. The blue line shows the slowest doubling time 
necessary to explain the raising part of the flare at C.L. of 
95%. The figure is reprinted from [6]. 


with a constant reveals a x 2 /N.d.o.f of 199/58 corre¬ 
sponding to a probability of 2.6 x 10 -17 . 

We use the rapidly rising part of the 1 st big flare 
(MJD 56244.0620-56244.0652) in order to compute 
the conservative, slowest doubling time, 7 d, which is 
still consistent with the MAGIC data. We fit the light 
curve with a set of exponential functions, each time 
assuming a given td value and computing the corre¬ 
sponding fit probability. We obtain that 4.9 min is 
the largest value of td, which can still marginally fit 
the data with probability > 5% (see the blue line in 
Fig. [2]). Note that the corresponding time scale in the 
frame of reference of IC 310 will be slightly shorter: 
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FIG. 3: MAGIC measurement of average spectral energy 
distributions of IC 310 during the flare (red) . For compar¬ 
ison we show the results from the high (blue, open squares) 
and low (black, open markers) states reported in [f| and 
the average results (gray triangles) reported in 0]. The 
dashed lines show power-law fits to the measured spectra, 
and the solid line with filled circles depicts the spectrum 
corrected for absorption in the extragalactic background 
light according to [12] ■ As a reference, the spectral power- 
law fit of the Crab Nebula observations from 4 is shown 
(gray, solid line). The figure is reprinted from [g]. 


FIG. 4: Significance map (color scale) of the Perseus 
cluster in gamma rays observed in the night of Novem¬ 
ber 12/13 th , 2012, with the MAGIC telescopes. The inset 
shows the radio jet image of IC 310 at 5.0 GHz obtained 
with the European VLBI Network (EVN) on October 29, 
2012. Contour lines (and associated to them color scale) 
increase logarithmically by factors of 2 starting at three 
times the noise level. The ratio of the angular resolution 
between MAGIC and the EVN is 1:580 000. The figure is 
reprinted from [g]. 



4.9/(l + z) min = 4.8 min. 

The observed spectrum can be described by a simple 
power law (see Fig. [3|: 

§ = (1) 

The flux normalization at 1 TeV obtained from the fit 
is f 0 = (17.7±0.9 stat ±2.1 syst ) x 10 _12 TeV _1 cm -2 s -1 
Even while the mean flux during the flaring night is 
4 — 30 larger than previous measurements, the spec¬ 
tral index, T = 1.90 ± 0.04 stat ± 0.15 syst is consistent 
with them within the statistical and systematic errors. 
No significant bend or cut-off is seen in the spectrum 
up to TeV energies. As part of the observation was 
carried out with a higher then usual offset angle from 
the camera center the systematic error on the flux nor¬ 
malization is slightly larger (12%) then reported in Q. 
The error of the energy scale is 15%. 


peak flux density of 77 mJy/beam and a lcr noise level 
of 0.027 mJy/beam. The restoring beam has a major 
and minor axis of 4.97 x 1.24 mas 2 with the major axis 
at a position angle of —8.5°. It contains a total flux 
density of Stot a i = 109 mJy, which we conservatively 
assume to be accurate to 10%. The dynamic range 
DR of the image, defined as the ratio of the peak flux 
density to tripled noise level in the image is ~ 950. 

The angle 6 of the radio jet to the line-of-sight can 
be determined from Doppler boosting arguments for a 
given jet speed /3 and spectral index a by considering 
the ratio R of the flux density in the jet and counter¬ 
jet: 


/1 + f3 cos 0 \ 
\ 1 — /3 cos 0 ) 


( 2 ) 


Following [nil ] we use the DR as an upper limit for 
the detection of a counter-jet. This gives us an upper 
limit of 0: 


B. EVN 


0 < arc cos 


(DR}'O-o0 _ i\ 

\DR x R 2 -“) + 1 ) ' 


(3) 


IC 310 has been observed with the EVN at 1.7, 5.0, 
8.4 and 22.2 GHz between 2012-10-21 and 2012-11- 
07. The data were amplitude and phase calibrated 
using standard procedures with the Astronomical Im¬ 
age Processing System (AIPS, 14]) and imaged and 
self-calibrated using DIFMAP 28f . 

In inset panel of Fig. 0] we present the image with 
the highest dynamic range obtained from the obser¬ 
vation at 5.0 GHz from 2012-10-29. The image has a 


Substituting DR in Eq. (2[ assuming a flat spectral 
index of a = 0 and /3 —>• 1, we obtain an upper limit 
for the angle between the jet and the line-of-sight of 
< 20 °. 

Additionally, the extension of the projected one¬ 
sided kpc radio iet of ~ 350 kpc measured at a wave¬ 
length of 49 cm [29j yields an estimate of a lower limit 
for the angle. De-projecting the jet using the upper 
limit quoted above would results in a lower limit of the 
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jet length of ~lMpc. Radio galaxies typically show 
jets extending up to 150kpc-300kpc |24|. The max¬ 
imal length of radio jets has been measured to be a 
few Mpcs which corresponds to an angle of ~ 5 — 10° 
in the case of IC310. Smaller angles would rapidly 
increase the de-projected length of the jet to values 
far above the maximum of the distribution of the jet 
lengths. 

III. INTERPRETATION 

GeV and TeV gamma-ray emission from blazars and 
radio galaxies is often explained in terms of shock-in- 
jet models. Charged particles are accelerated in an 
active region moving along the jet. Causality con¬ 
dition provides that the variability time scale of the 
observed emission can be used to constrain the size of 
the emission region. 

A conservative estimate of the shortest variability 
time scale in the frame of reference of IC310 yields 
Ai/(1 + z) = 4.8 min. Using the best mass estimate 
of IC310 black hole this measurement corresponds to 
20% of the light travel time across the event hori¬ 
zon. Even allowing for the factor 3 uncertainty in 
the mass the fraction, 60%, is still below one. The 
ultrafast variability casts a shadow of doubt on the 
current shock-in-jet paradigm. The moving shock 
plasma leads to a shortening of the observed variabil¬ 
ity time scale At compared with the variability time 
scale At' in a frame comoving with the shock given by 
At = (1 + z)5~ 1 At'. This effect is often used to ex¬ 
plain ultrafast variability from blazars Hi in which 
5 can be nearly arbitrarily large providing that the 
jet moving with large Lorentz factor is observed at a 
very small angle. In the case of IC 310 however the es¬ 
timation of the observation angle 10° — 20° obtained 
from the radio observations constrain the maximum 
Doppler factor to be < 6. 

All of these attempts to explain the sub-horizon 
scale variability with relativistic projection effects 
alone encounter a fundamental problem [20|. If the 
perturbations giving rise to the blazar variability are 
injected at the jet base, the time scale of the flux vari¬ 
ations in the frame comoving with the jet is affected 
by time dilation with Lorentz factor Tj. In blazars 
where <5 ~ Tj, the Lorentz factor cancels out, and the 
observed variability time scale is ultimately bounded 
below by AteH- 

Additionally, a very high value of the Doppler 
factor is required to avoid the absorption of the 
TeV gamma rays due to interactions with low- 
energy synchrotron photons. Such synchrotron 
photons are inevitably produced together with 
the gamma rays in the shock-in-jet scenario. The 
optical depth to pair creation by the gamma 
rays can be approximated by r 77 (10 TeV) ~ 
300 (5/ 4) -6 (At/1 min) -1 (L syn /10 42 erg s _1 ). 


Adopting, conservatively, a non-thermal infrared lu¬ 
minosity of ~ 1% of the gamma-ray luminosity during 
the flare, the emission region would be transparent to 
the emission of 10 TeV gamma rays only if 5 > 10. 

In summary, trying to interpret the IC 310 flare in 
the framework of the shock-in-jet model meets diffi¬ 
culties. Alternative models can involve stars falling 
into the jet II, mini-jet structures within the jets 
pj| or magnetospheric models Ed, ns m HI- In the 
case of IC 310 star-in-jet model cannot provide suf¬ 
ficient luminosity to explain the TeV flare [f|. Also 
jets-in-jet models suffer from rapidly dropping lumi¬ 
nosity at larger observation angles |6j. Moreover the 
magnetic reconnection which can led to production of 
such mini-jets is expected to occur in the main jet 
rather at larger distances from the black hole. 

In magnetospheric models, particle acceleration is 
assumed occur in electric fields parallel to the mag¬ 
netic fields. This mechanism is common to the 
particle-starved magnetospheres of pulsars, but it 
could also operate in the magnetospheres anchored to 
the ergospheres of accreting black holes (see Fig. [5]). 
Electric fields can exist in vacuum gaps when the 
density of charge carriers is too low to cause their 
shortcut, i.e. below the so-called Goldreich-Julian 
charge density. Electron-positron pairs in excess of 
the Goldreich-Julian charge density can be produced 
thermally by photon-photon collisions in a hot accre¬ 
tion torus or corona surrounding the black hole. It 
has also been suggested, that particles can be injected 
by the reconnection of twisted magnetic loops in the 
accretion flow [22j . A depletion of charges from ther¬ 
mal pair production is expected to happen when the 
accretion rate becomes very low. In this late phase 
of their accretion history, supermassive black holes 
are expected to have spun up to maximal rotation. 
Black holes can sustain a Poynting flux jet by virtue of 
the Blandford-Znajek mechanism [ll|. Jet collimation 
takes place rather far away from the black hole, i.e. 
at the scale of the light cylinder beyond ~ 10r g . Gaps 
could be located at various angles with the jet axis 
corresponding to the polar and outer gaps in pulsar 
magnetospheres leading to fan beams at rather large 
angles with the jet axis. As the gap height and seed 
particle content depend sensitively on plasma turbu¬ 
lence and accretion rate, the gap emission is expected 
to be highly variable. For an accretion rate of ~ 10 -4 
of the Eddington accretion rate and maximal black 
hole rotation, the gap height in IC 310 is expected to 
be h ~ 0.2r g da which is in line with the variability 
times seen in the observations. Depending on the elec¬ 
tron temperature and geometry of the radiatively inef¬ 
ficient accretion flow, its thermal cyclotron luminosity 
can be low enough to warrant the absence of pair cre¬ 
ation attenuation in the spectrum of gamma rays. In 
this picture, the intermittent variability witnessed in 
IC 310 is due to a runaway effect. As particles accel¬ 
erate to ultrahigh energies, electromagnetic cascades 
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FIG. 5: Scenario for the magnetospheric origin of the 
gamma-rays: A maximally rotating black hole with event 
horizon r g (black sphere) accretes plasma from the cen¬ 
ter of the galaxy IC 310. In the apple-shaped ergosphere 
(blue) extending to 2r g in the equatorial plane, Poynting 
flux is generated by the frame-dragging effect. The rota¬ 
tion of the black hole induces a charge-separated magne¬ 
tosphere (red) with polar vacuum gap regions (yellow). In 
the gaps, the electric field of the magnetosphere has a com¬ 
ponent parallel to the magnetic field accelerating particles 
to ultra-relativistic energies. Inverse-Compton scattering 
and intense pair production due to interactions with low- 
energy thermal photons from the plasma accreted by the 
black hole leads to the observed gamma rays. The figure 
is reprinted from [gj. 

develop multiplying the number of charge carriers un¬ 
til their current shortcuts the gap. The excess parti¬ 
cles are then swept away with the jet flow, until the 
gap reappears. 

IV. CONCLUSIONS 

Radio galaxies and blazars with very low accretion 
rates allow us to obtain a glimpse of the jet forma¬ 
tion process near supermassive black holes. Observa¬ 


tions of IC 310 performed with the MAGIC telescopes 
showed variability with time scale below 5 min, shorter 
than the light crossing time of the event horizon of its 
black hole. The commonly used in AGNs shock-in- 
jet models have troubles to explain such emission. A 
plausible explanation involves emission from vacuum 
gaps in the magnetosphere of IC310. Interestingly, 
such explanation invite to explore analogies with pul¬ 
sars where particle acceleration takes place in two 
stages. In the first stage, particle acceleration occurs 
in the gaps of a charge-separated magnetosphere an¬ 
chored in the ergosphere of a rotating black hole, and 
in a second stage at shock waves in the force-free wind 
beyond the outer light cylinder. 
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